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Abstract 



l/-^ ■ Wavy film flow of incompressible Newtonian fluid down an inclined plane is 

o. 

^^ ■ considered. The question is posed as to the parametric conditions under which 



the description of evolution can be approximately reduced for all time to a 
single evolution equation for the film thickness. An unconventional perturba- 
tion approach yields the most general evolution equation and least restrictive 



^ ■ conditions on its validity. The advantages of this equation for analytical and 

numerical studies of 3D waves in inclined films are pointed out. 
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I. INTRODUCTION 

Thin liquid layers ("films") flowing along solid surfaces — such as inclined planes or ver- 
tical cylinders — occur in both natural and man-made environments. Industrial applications 
of fllm flows started as long ago as the 1800s and have been growing in their scope and 
importance ever since (see e.g. Refs. |I|J3). 

Accordingly, the studies of fllm flows, in particular those down an inclined plane, have 
a considerable history (see e.g. Refs. |T],^. However, the dynamics of nonlinear waves (typ- 
ically present in fllm flows) is still far from being satisfactorily understood (see Refs. ^J^ 
for the most recent progress reviews). The nonlinear Navier-Stokes (NS) partial differential 
equations (PDEs) of this "Kapitza problem" couple together several flelds (pressure and the 
components of velocity), each a function of time and three spatial coordinates. Furthermore, 
the boundary conditions (BCs) of the problem involve a free boundary (the free surface of 
the fllm) whose PDE itself is coupled to the NS equations. The full spatially 3-dimensional 
(3D) problem is too hard to simulate even with the most powerful modern computers. 

Even simpler, 2D computations of wavy fllms have been undertaken only under the 
simplifying assumptions of short spatial intervals (e.g. Refs. ^^ and/or time-independence 
(e.g. Ref. ^. However, 2D flows are frequently unstable to 3D disturbances, and three- 
dimensionality can be important for many inclined fllms (see e.g. recent experiments in Ref. 

Therefore, naturally, one looks for more manageable approximate descriptions of wavy- 
fllm evolution. Such simplifled theories are possible in certain domains of the space of pa- 
rameters (for which, typically, the slopes of the surface waves are small). Of course, the 
greater simpliflcation has to be paid for by a more limited applicability of the theory. The 
greatest simpliflcation is achieved when the problem reduces to a single evolution PDE ap- 
proximating the thickness of the fllm. [From the numerical-simulation point of view, the 
most important simpliflcation here is the reduction in the number of independent spatial 
variables, and even a lower- dimensional system of equations would have been not much more 



difficult than a single evolution equation (EE). However, no such system has ever been de- 
rived consistently for an inclined single layer film, the subject of our consideration in this 
paper. The same is true for the nonlocal, i.e. integro-differential, film-thickness equations. 
(The only known examplet^ of such a nonlocal EE, as far as film fiows are concerned, is 
an EE for a core-annular fiow, and not for an inclined film.) In any case, in this paper, the 
nonlocal equations are altogether excluded from the consideration; thus, when we say "the 
most general EE" it should be understood as "the most general local EE" , etc.] Although 
less drastic simplifications than a single EE, and therefore having a larger parametric range 
of validity, are known (see e.g. Ref. ^, so far fully- dimensional simulations for sufficiently 
extended spatial domains have been carried outall^ only for the theories hinged on a single 
evolution equation (EE). Such single-EE theories of inclined-film fiows are the subject of the 
present paper. 

Evolution equations for film thickness have been known since the pioneering work of 
Benney.0 The conventional perturbation approach to their derivation (e.g. Refs. ^M3 
used a small (longwave) parameter, say e. In particular, each of the "global" ("internal", 
"basic" ) parameters specifying the problem (the parameters appearing in the dimensionless 
NS equations and the free-surface BCs) must be ascribed, in such a single-parameter (SP) 
technique, a certain power of e as its order of magnitude (Om)- Therefore, an artificial depen- 
dence is forced on the — intrinsically independent — parameters. This unnecessarily restricts 
the domain of justified validity of the resulting EE. For example, for the vertical film, there 
are just two independent parameters: the "Reynolds number" R and the "Weber number" 
W. If R is of the order of magnitude of (~) e" and W ~ e^, then W ~ R^^"-. The set of 
points W = R^f'^ is just a ID curve in the 2D space (i?, W), while the complete domain for 
which the EE is valid is likely to have the same dimensionality as the parameter space (i?, 
W) itself, that is it should be a 2D domain. 

Furthermore, each time the powers assigned to the parameters are altered, it is in prin- 
ciple necessary to again go through the entire procedure or the SP derivation, and one can 
arrive at a different EE as a result. For example. Topper and Kawaharap3 considered two 



cases of an inclined-film flow (such a flow is specified by three global parameters: in addition 
to R and W, the inclination angle 6 can be independently varied). In one case, they required 
the angle (p (= 7r/2 — 6) of the film plane with the vertical to be small, ~ e, and stipulated 
R r^ e^ and W ~ e~^ (in our definitions of R and W; see section y below). As a result, they 
obtained an EE containing both dissipative and dispersive terms (and with all coefficients in 
the EE being ~ 1). However, their second case, for which they chose cot^ ~ 1, /? ~ 1 and 
W ~ e~^, resulted in an equation with no dispersive terms. (In the SP framework, it is only 
formally that the latter equation can be obtained from the former by omitting its dispersive 
term, and the only way to really "justify" the nondispersive equation is to repeat the entire 
derivation procedure starting all the way back from before NS equations.) 

If a given inclined-film system is not close to any of these two parameter-space curves, 
the theoryllB is invalid. Logically, there are three possibilities: (i) the flow evolution cannot 
be (approximately) reduced to a single EE for all time; (ii) such a reduction is possible but 
the resulting EE is different from each of those obtained in Ref.E^; or (iii) the EE coincides 
with one of their two EEs, despite the different parameter curve. We are naturally led to the 
following questions: (i) under what (parametric) conditions is an approximate description 
of the film flow possible (for all time) which can be reduced to a single EE? (ii) How can 
the set of all such EEs be characterized? These are the questions we pose and attempt to 
answer, for inclined-film flow, in this paper. 

We note that one must distinguish between the all-time validity of an EE and the limited 
in time validity. This issue arises because in the SP approach, the fixed powers of the 
small parameter are stipulated not only for the global parameters of the system, which do 
not depend on time, but also for the characteristic time- and length-scales. However, these 
characteristic scales can change with time as the long — dissipative — system proceeds to the 
attractor; they are instantaneous parameters. So, the assumption of fixed scales is incorrect 
after a limited time has elapsed, and therefore, the SP derivation is invalid. 

The rest of the paper is organized as follows. In section II, we formulate the full NS 
problem. In section III, we introduce an iterative perturbation procedure. It starts with the 



well-known (although typically unstable), waveless, "Nusselt" solution of the NS problem. 
The only principle necessary in deciding the iteration steps is the requirement that in the end, 
a single (and valid for all time) EE should be arrived at, with minimal simplification of exact 
equations. No dependencies are imposed on the internal parameters: the validity conditions 
(VCs) we obtain require that several quantities, which are certain products of powers of 
the internal parameters, be small — independently of one another. [So, several independent 
(small) parameters emerge in the derivation. Thus, the iterative procedure we introduce 
here is a variation of the "multiparameter" perturbation (MP) approach developed in Refs. 
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In section IV, we arrive at the most general evolution equation (GEE) valid (provided 
certain restrictions on parameters are satisfied) for all time: Any all-time-valid EE derivable 
with a conventional single-parameter (SP) approach necessarily coincides with one of just 
a few prototype equations which are certain truncations of the GEE. The GEE also has 
the least restrictive domain of validity: It is easy to obtain the domain of validity for each 
simplified EE, and it is a subdomain of the all-embracing (corresponding to the GEE) 
validity domain. Outside of the latter parametric domain, there is no single EE which could 
approximate the film evolution for all time. That this is the case is argued in section V 
by analyzing the structure of possible correction terms of the EE through all orders of the 
iteration procedure. (In particular, as was shown in Ref. |2l| and Ref. ^2| — where a different 
derivation of the same EE was sketched — the amplitudes of inclined-film waves which can 
be described (for all time) by a single EE are necessarily small.) The paper is summarized 
in the last section. Some of the more technical considerations are relegated to Appendices. 

Simulation results for different versions of the inclined-film flow EEs we obtain here will 
be given in a sequel to this paper. [For some of those results (in particular, those showing 
good agreement with experiments^) , see Ref. ^.] 



II. THE EXACT NAVIER-STOKES PROBLEM 

We consider a layer of an incompressible Newtonian liquid flowing down an inclined plane 
under the action of gravity. Our (cartesian) coordinates are as follows: the x axis is normal 
to the plane and directed into the fllm; the y axis is in the spanwise direction; and the z axis 
is directed streamwise (the overbar here and below indicates a dimensional quantity). The 
corresponding components of velocity are u, v, and w. We denote by p the pressure fleld in 
the fllm; the pressure of the ambient passive gas is neglected for simplicity. 

The system is determined by the following independent (dimensional) parameters: the 
average thickness of the fllm /iq; the hquid density p, viscosity /I, and surface tension a; 
gravity acceleration 'g] and the angle of the plane with the horizontal 6. 

There is a well-known time- independent, Nusselt's solution of the NS problem for the 
inclined fllm. The thickness of the Nusselt fllm is constant (hence, Nusselt's flow is also re- 
ferred to as a "flat-fllm" solution). The only nonzero component of velocity is the streamwise 
one. It only changes across the fllm, starting from the zero value at the solid plane. The 
free-surface value U of the Nusselt velocity isU = ghl sin /{2u) (where u = ]!/ pis the kine- 
matic viscosity). We nondimensionalize all quantities with units of measurement based on 
p, ho, and f7 (we have, e.g., the following units: h^ for all coordinates; U for velocities; ho/U 
for time t; pU'^ for pressure; pU'^ho for surface tension; etc.). We will see that exactly three 
independent basic parameters (BPs) appear in the dimensionless equations and boundary 
conditions of the problem; one can choose e.g. the inclination angle 9, the Reynolds number 
R = hoU/u [= ghlsin6/{2i?'^)], and the Weber number W = aR/2 [= a/{pghlsm6)], as 
such BPs. 

We can write the Navier-Stokes momentum equations in coordinates moving relative to 
the solid plane with a constant velocity V in the z-direction (i.e. introducing z = z — Vt 



and omitting the tilde) in the following dimensionless form (see e.g. Ref. pSj) : 

Ut + UUx + VUy + WUz — VUz = 



2 1 

PX- ^COte + —{U^^+Uyy+U^z), (1) 



Vt + UVx + VVy + WVz — VVz 



Py + ^ (^xx + Vyy + Vzz), (2) 



Wt + UWx + VWy + WWz — Vwz = 

2 1 

-PZ + -f^ + -f^{WxX+Wyy+Wzz)- (3) 

K K 

(The subscripts x, y, z, and t here and below denote the corresponding partial derivatives. 
We will see below that it is appropriate to choose V = 2, the common phase velocity of all 
infinitesimally weak waves.) The continuity equation is 

Ux + Vy + Wz = 0. (4) 

The BCs are as follows. The no-slip conditions at the solid plane are 

u = V = w = (x = 0). (5) 

The tangential-stress balance conditions at the free surface x = h{y,z,t), the local film 
thickness, are 

{Vx + Uy){l - hy) + 2{Ux - Vy)hy - {v ^ + W y) k ^ 

- {uz + Wx)hyhz = {x = h) (6) 
and 

{uz + Wx){l - hi) + 2{ux - Wz)hz - {vz + Wy)hy 

- {uy + Vx)hyhz = {x = h). (7) 



The normal-stress balance condition is 



2 

-p{i + hl + hlf^ ~ ;r h=^ + ^yK + ^^^ 



~{Uy + Va:)hy + {Vz + Wy)hyh 

1/2 



-(u, + u;,)/i,](l + /ij + /i2) =a 



/ij^j^(l + h^ 



+h^^{l + hy ) - 2hyh^h^y'^ {x = h). (8) 

Finally, the kinematic condition at the free surface is 

ht + vhy + whz — Vhz = u {x = h). (9) 

III. ITERATIVE PERTURBATION PROCEDURE 

A. Minimal requirement of derivability 

As was motivated in the Introduction, we are interested in the question of (approximate) 
reducibility of the above complicated description of inclined-film dynamics to a single EE. 
It appears that an iterative perturbation approach is appropriate for the general analysis of 
the problem. 

Looking at the previously known, conventional (single-parameter expansion) derivations 
of EEs (valid each for its own particular curve in the parameter space; see the Introduction), 
it is clear that each of them in effect discards some terms of the NS momentum equations so 
that each of those essentially becomes an ordinary differential equation (ODE) in x, linear 
and with constant coefficients, that is easy to solve (with similar simplifications in BCs). 
When these solutions (for velocities in terms of film thickness) are substituted into the 
kinematic condition (^, the EE for thickness h (or, equivalently, for the thickness deviation 
1] = h — 1) ensues. However, after all the quantities have been expanded in power of e (see 

8 



Introduction), one has no control over which NS terms to omit: this is simply dictated by 
the expansion scheme of the SP approach. Thus, sometimes "harmless" terms are discarded: 
even if they were retained, one still would be able to solve for velocities, arriving, as a result, 
at a clearly more general EE. 

Accordingly, our main idea in this paper is to look for a derivation in which the single 
postulated requirement would be that a maximally general EE for film thickness be arrived 
at in the end; so only those terms of the exact NS equations will be discarded which are 
clearly in the way of obtaining linear ODEs for velocities and pressure [we call this the 
"minimal requirement of derivability" (MRD)]. In this way we obtain approximate solutions 
for the deviations of exact solutions from the "seed" Nusselt's fields [e.g., the approximation 
Wq to Wq = w — Wj\f, where the Nusselt Wn is known: see (|12D below]. For some parameter 
values, this immediately leads to an EE whose approximation of exact evolution is good 
for all time (here, as was mentioned in the Introduction, we are only interested in such 
all-time-valid EEs; see e.g. Ref. ^ for a further discussion of their difference from the EEs 
whose validity is limited in time). In other cases, however, the procedure must be repeated, 
with the refined solutions, such as w^ + Wq, playing the seed role that was played by the 
Nusselt solutions on the original stage (so that at the second iteration stage one determines 
the approximation wi to wi ^ w — wn — wq, etc.). Thus, ours is an iterative perturbation 
approach — which is known even in general to be an alternative to the expansion method 



(see e.g. Ref. |2|). 

Estimating all members of NS problem equations in terms of parameters, the requirement 
that the discarded members be much smaller than those retained leads to the parametric 
conditions for our derivation to be valid (see Appendix ^. Thus, the method yields (i) the 
evolution equation for film thickness; (ii) explicit expressions for velocity components and 
pressure in terms of the film thickness; and (iii) the parametric validity conditions of the 
theory. 

It is known (see Ref. ^ and references therein) that no single-EE description can exist 
globally in time for those parametric regimes of inclined-film flow which lead to the eventual 



amplitude of surface waves being "large" — comparable to the average film thickness. But in 
the present communication, as was mentioned above, we are interested exactly in the large- 
time behavior, when the system is already close to the attractor, and we want a single-EE 
description of the wavy film dynamics. Therefore, in addition to the MRD, we can use from 
the very beginning — in order to simplify our derivation — the requirement that the amplitude 
A of the film thickness deviation, A{t) = max \ri\ , be small (for all time): with 

h=l + r], (10) 

we have 

mixx\r]{y,z,t)\ = A<^1. (11) 

[Note that A can depend upon time; in such cases, we say that the parameter is a local 
parameter, in contrast to the time-independent, "global" BPs 6, R, and W.] However, unlike 
the conventional derivations, we do not have to postulate that the characteristic lengthscales 
in the film plane are large (the longwave assumption); rather, this will follow from the MRD. 

B. Nusselt's solution 

The above-mentioned Nusselt solution (of the NS problem) which is steady and uniform 
along the film (i.e. in the streamwise and spanwise directions) is as follows. The dimensionless 
Nusselt streamwise velocity wn is 

wn{x) = 2x — x^ . (12) 

This clearly satisfies the z-NS equation 

wnxx = -2 (13) 

with the boundary conditions 

WNx = (x = l) (14) 

10 



and 

wn = {x = 0). (15) 

Similarly, the Nusselt pressure 

Pn = ^icote){l-x) (16) 

is the solution of the x-NS equation 

2 

Pnx = --f^ cot 9 (17) 

K 

with the boundary condition 

Pn = (x = l). (18) 

Finally, the Nusselt normal and spanwise velocities are 

Mat = and vjy = 0. (19) 

At sufficiently large Reynolds number, the destabilizing effect of inertia overcomes the 
stabilizing influence of gravity so that the Nusselt solution loses its stability. (An analysis 
later in this section leads to R > (5/4) cot^, the well-known criterion for instability.) As a 
result, the film is not uniform any more and also changes with time. 

C. First iteration step 

We represent the exact velocities and pressure in the form of sums of Nusselt 's solutions 
and the deviations from those: 

W = Wn + Wq, U = Un + Mo, 

V = vn + Vo, p = p^+ po, (20) 
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where tildes indicate the deviations from the Nusselt solutions. First, we consider the z- 
momentum NS equation. We substitute the expressions (^) for velocities and pressure into 
Eq. (^. The resulting (exact) equation can be written in the form 

W^xx = RpOz - V^Wo + RWot + RUoiwN + Wo)x 

+ Rvowoy + R{wn + wo- V)woz, (21) 

where V^ = d'^/dy'^ + d"^ /dz^ . In accordance with the MRD, as was discussed above, in 
order to obtain a solvable ODE in x, we have to discard all the terms on the right-hand-side 
(RHS) of (0), since every one of those contains unknown quantities. This yields a simplified 
equation for the velocity Wq, whose solution we denote by Wq: 

woxx = 0. (22) 

We call Wi the error in the approximation of Wq by Wq. 

Wq = Wq + Wi. (23) 

It is clear [see (|TD|) and (|T1|)] that the characteristic x-lengthscale X is of Om of 1 (X ~ 1), 
i.e. (since d/dx ~ 1/X), 

We denote by Y and Z the characteristic lengthscales in the spanwise and streamwise di- 
rections respectively, so that 

d I 

25) 

dy Y ^ ' 

and 

d I 

IT. ~ r (^*" 

Similarly, denoting by T the characteristic timescale, we have 

12 



- - -• (27) 

dt T ^ ^ 

Also, clearly, d^dx^ ~ 1, d^dz^ ~ 1/Z^, d/dy^ ~ 1/Y^, etc. 

The characteristic magnitude of each of the neglected terms on the RHS of Eq. (|2l|) — 
estimated by replacing wq with wo, etc. — should be smaller than that of the term on the 
LHS, Wqxx'- Wq/Z"^ ^ Wq/X"^ and Wq/Y"^ <^ wq/ X"^. Hence, we have the following restrictions 
on the (instantaneous) lengthscales for the validity of our theory 

^ « 1 (28) 

and 

i « 1- (29) 

Thus, as a consequence of our derivability requirement, we have obtained the longwave 
conditions — which are rather postulated in the conventional "lubrication" or "long-wave" 
derivations. 

Similarly, using (pGf ) and (pT]), the conditions that the advective inertial term contain- 
ing wq, RVwqz, and the time- derivative Rwot be smaller than Wq^x lead to the following 
requirements on the parameters: 

XL , , 

- « 1 (30) 

and 

§ « 1- (31) 

Henceforth, we confine our consideration to the film configurations which satisfy the condi- 
tions (p8D-(pID. Note that since the local parameters A, Z, Y, and T can change with time, 
they may cease to satisfy the validity conditions [such as ([TT| ) and (p^)-(|3T|)] at a certain 



stage of the evolution. In such cases, clearly, the single- EE description would be valid only 
for a limited time that these (and some other, additional conditions appearing in Appendix 

13 



^ still hold. Later, we will determine domains in the space of global parameters for which 
such a violation of VCs never happens — so that the single-EE approximation is valid globally 
rather than merely locally in time. 

The boundary condition on wq, at x = /i, given by the tangential-stress balance equation 
(0), is written in the form 

Wox = -wnx - Uo^ + [2(wo2 - uox)ri^ + {vq^ + WOy)Vv 

+ {uoy + vox)VyVzKl - V^y^ [x = h) . (32) 

According to the MRD, we have to drop all the terms containing the unknown quantities 
[those denoted by letters with tilde] on the RHS. Also, using the smallness of the surface 
deviation (pTf) , we will everywhere transfer the boundary conditions from the true boundary 
X = hio a. convenient "boundary" x = 1 by expanding all quantities in Taylor series around 



X = 1 (cf. Ref. ^^, such as 

wox{x = h)= wox{x = 1) + woxx{x = 1)77 H . (33) 

Noting that wn{x = 1) = 0, we have the simplified boundary condition (for the approxima- 
tion wq of the exact quantity wq): 

wox = -wnxxV = 2t] (x = 1). (34) 

Finally, the no-slip condition requires 

Wo = {x = 0). (35) 

The solution of Eq. (^) with the boundary conditions ( p^ and ( pSf ) is clearly 

Wo = '2r]x. (36) 

Note that from (|g) and (|3§) 

wnx + Wox = {x = h), (37) 

14 



which we will use below. We also observe that 

u7o ~ ^ (< wjv ~ 1). (38) 

Next, consider the incompressibility equation [see Eq. @] in the form 

UOx = -VOy - Woz - Wi^- (39) 

(Note that wn does not make any contribution to this equation, since wjy does not depend 
on z.) Dropping the terms with unknowns vq and Wi on the RHS, we obtain an equation for 
the approximation Uq: 

Uox = -Wq:, = -2X7]^. (40) 

The no-shp BC (||) requires 

Mo = {x = 0), (41) 

and one readily obtains the solution 

Uq = -X^T]^. (42) 

We note that 

uo ~ - < 1. (43) 

We denote the deviation of the exact solution uq from the approximation uq by Mi, so that 

uq = uq + ui. (44) 

The (x-momentum) NS equation for the deviation of pressure from the Nusselt solution 
Pn is 

I, . I ,. ^ ^ . 

POx = -^{Uoxx + UOyy + Uq^z) + -^{Ulxx + Uiyy + Mi^J 

JrC JrC 

-{Uq + Ui)t - {Uq + Ui){Uq + Ui)x - Vq{Uq + Ui)y 

15 



- {wn + wo + wi-V){uo + ui)^ (45) 

[where we have used Eq. (0)]. In this equation, we have to drop the terms containing the 
unknowns Vq, Wi, or Ui. Next, it is easy to show that the Uq^x term [which is ~ Uq due to 
(1^)] on the RHS is much larger than the other terms. Namely, Uoyy ~ Uq/Y"^ <^ Mq ~ Uq^x 
due to (1^) and uq^z ~ Uq/Z"^ ^ Mq ~ Uq^x due to (|29|). Similarly, Ruot ~ {R/T)uq <^ 
Uq ~ Uqxx by making use of (0) and R{wj\f + Wq — V)uqz ~ Ruqz ~ {R/Z)uq ^ Mq ~ Moxi 
due to (|30D . Finally, uqUqx -C mqx ~ moxx, since mq <^ 1 [see Eq. (|43|) ]. Thus the simplified 
equation (for the approximation pq) is 

1 2 

Pox = "H^oxx = -^Vz- (46) 



o 
PO = -PN + -^ [[{Uo + Ul)x + VOyVy + {wq + Wl)zril - [{uq + Ui)y + Vox]r]y + [vqz + {wq + Wi)y]r]yr]z 



The BC at X = h [see the normal-stress balance condition, Eq. (^] is 

2 
R 

-[{uq + ui)z + Wix]Vz} (1 + ^J + vly^ - ^ [Vyy (l + ^^) + Vzz (l + r^j) - 2r]yr]z%z 

X (l + V'y + Vl)-'^' ix = h) (47) 



[where we have used (|3^) to eliminate wnx and ti^oa;] • We have to drop those terms containing 
Vq, Wi, or Ml. In addition, a number of known terms are estimated to be smaller than the 
term with Uq^- Namely, by using the estimates of Wq and Uq, (^) and (p3D, WQyrjyTjz ~ 
{A/Z){A/Yf < A/Z ~ Mo ~ Mox- Similarly, wo.r^^ ~ {AIZ){A^ /Z'^) < mox- Also, mo?;^/^ ~ 
Uov4/F^ ^ Uox and UqzVz ~ Uq^/^^ ^ ^ox- Finally, one can see that each of the other 
known nonlinear terms can be neglected since r/^ ~ A? jZ"^ ^ 1, r\^ ~ A? jY''^ <^ 1, and 
VyVz^yz ~ A^ /{Y'^Z'^) <^ 1. We also note that the entire Taylor expansion for p^i^x = l + rf) 
consists of just one term, 

PN{h)=PNx{l)ri- (48) 

With the known expressions for p^, wq, and uq (see above), and using Taylor series about 
X = 1, truncated to the first nonzero term, for all the quantities, one obtains the BC 
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2 ^2 



(cot e)r] - -T], - aV^r] (x = 1). (49) 

K K 



The solution of Eq. (46) with the boundary condition (^9[) is 

Rpo = 2(cot 9)ri - 27^^(1 + x) - 2WV'^t]. (50) 

We observe that 

f ^ I W W\ ^ ,,^, 

i?Po~max I cot 6', — , — , ttj ) ^- (51) 

As usual, we call pi the difference between pQ and its approximation pq: 

Po=Po+Pi- (52) 

Finally, consider the ?/-NS equation, 

't'Oxi = R{PO +Pl)y - V^{)o + Rvot + R{uq + Ui)vo^ 

+ RvQVQy + R{wn + wq + wi-V)vqz- (53) 

Dropping the unknown terms on the RHS, we obtain the simplified equation: 

"VQxx = RPOy 

= 2(cot 9)r]y - 2(1 + x)r]y, - 2WV\, (54) 

where we have used the expression (^) for po iii terms of r] to get the RHS in the explicit 
form. 

The BC on vq, at x = h, is [see Eq. (|)] 

VQx = -(mo + Ul)y + {2[t^0y 
-(UO + Ui)x]Vy + [{Uo + Ui)z + Wlx]VyVz 
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+ [voz + (wo + wi)y]ri^} (l - 7]fj {x = h) 



(55) 

where we have again used ( p7| ) to ehminate wnx + wqx- Omitting the terms with tildes and 
performing the by now famiUar estimates of the known terms on the RHS, we see that the 
term uoy is larger than any other term, so that the simplified BC at x = 1 is 



Vqx — ~Uoy — r]zy 



(x = l), 



(56) 



where the expression (^21) for uq has been used. Finally, the no-slip condition 



vo = Q (x = 0) 



(57) 



is to be satisfied. The solution of Eq. (|5^ ) with the BCs ([56|) and (|^) — as can be readily 
checked by direct substitution — is 



^^0 



(cot Oriy - WV'^riy) (x^ - 2x) 



X 



^zy\— +X -Ax 



This yields the following estimate for vq: 

"cot^ W W 1 
Vq ~ max 



y ' ys' y^2' yZ 



A. 



As usual, 



(58) 



(59) 



Vq = t-o + ^l- 



(60) 



Substituting ( |60l) into (^), we observe that in obtaining the simplified equation (|40| ) for 
Uqx, we have in fact discarded VQy. This requires 



^0 ^ 



(61) 



(which will be used below). By using here the estimates (P3| ) for uq and (^) for vq, we arrive 
at the VCs max [Zcot9/Y^, WZ/Y\ W/(Y^Z)] < 1. These are a subset of the complete 
set of (instantaneous) VCs obtained in Appendix ^ [see Eq. (|A10|) ]. It is straightforward to 



verify that each discarded term in every step of our procedure is small as a consequence of 
those VCs. 

The kinematic condition [Eq. (^] at x = h becomes 

Vt + {vo + Vi)riy + {wn + W0 + W1- V)r]^ 

= uo + ui. (62) 

Dropping the unknown terms with tilde and the smaller terms voriy ~ [vq/Y)!] <^ uq [see 
Eq. (|6lD] and woriz ~ A^/Z <^ uq, we have 

Vto + (^Tv - V)r]z = Uq, (63) 

where we have introduced the "fast" time to; such that dt = dt^ +dt^ ~ dt^. Using the Taylor 
expansions for wn and -uq, "we obtain, at x = 1, 

r/,, + (2 - V)r^^ = 0. (64) 

Choosing l^ = 2, we can eliminate the fast-time undulations (which are clearly due to the 
uniform translation of the wave, with no change in its shape): 

Vto = 0. (65) 

Thus, the leading approximation determines the velocity of a reference frame in which film 
thickness does not change on the fast time scale. However, it will change with the slower 
time ti. In order to obtain this slower-time evolution of the film thickness, one needs to 
consider the next approximation for the velocities and pressure. From now on, we fix 



V = 2. (66) 



Then, in view of (p5|), dt = dt^. 
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D. Second iteration 

We now proceed to consider the "corrections" Wi, Ui, pi, and Vi for the velocities and 
pressure. By substituting 

Wo = Wo + Wi, Uo = Uo + Ml, 

Vo = vo + Vi, Pq=pq+ P^ (67) 

into the z-NS equation (plD, and taking into account wqxx = (p2D, we have the exact 
equation 

Wlxx = -Woxx + RPOz - V^Wo + RWot 
+RUo{wn + Wo)x + RVoWoy + R{wn + Wq - 2)woz 

+ [terms containing Wi, Ui, Vi, or pi]. (68) 

Performing our standard simphfication procedure, i.e. discarding the unknown terms (con- 
taining tilde) and small terms, and also taking into account that wq <^ wn [see Eq. (pSDJ 
and Voy <^ Uq [see Eq. (0)], the simplified equation for the approximation wi is 

Wixx = RpOz - V^Wo + RuqWnx + RwnWqz 

-2Rwo, + Rwot = 2(cot Or], - WV%) - 2{x + l)r],, 

- 2xV^r] + {2x^ - Ax)R7], + 2xRrit (69) 

where we have used the known expressions (in terms of surface deviation rj; see the preceding 
section) for the first-iteration approximations Pq, Wq, and Uq. (For analogous equations of 
the general, nth iteration step, see Appendix y.) 
The BC on Wi, at a; = h, comes from (|^): 
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+ {uQy + VQ^)riyr]z 
+ [terms containing wi, -Ui, -Ui, or pi] (x = /i) (70) 



where we have taken into account (pTl) . Continuing to use the simphfication procedure es- 



tabhshed in the previous section, we arrive at the boundary condition 

wix = -uoz = Vzz (a; = !)• (71) 

All other terms on the RHS of Eq. ( [fOD are readily estimated in our usual way to be smaller 
than uqz [Eq. ( |6lD is useful in estimating terms containing v]. The no-slip condition is 

wi = (x = 0). (72) 

The solution of the problem (|69D, ([7TD, and (^) is 

wi = (^cot % - WV^r]z) (x^ - 2x) 

/x^ 2 3 4 \ / 2 . 2\ 

+ I — - -x + -x I i?r7^ + ( 5x - -X -X j rjzz 

+ (a^-yj ^s/j/+ (y -2;j /?r/i, (73) 

which can be verified by direct substitution into the problem equations. Note that all the 
terms of Wi are estimated to be quadratic in the local parameters ( |A10|) ; we will see that, 
in general (see Appendix |CD, i(;„ is of the power (ra + 1), and similarly for m„, t>„, and pn- 

Taking into account uqx = —w^z [see Eq. (^)], the incompressibility condition yields the 
equation for ui, 

Ulx = -Wiz - W2z - Voy - Viy, (74) 
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where we have expressed wi as wi = wi + W2. Dropping the unknown terms, we obtain an 
equation for the approximation Ui'. 

Ulx = -Wiz - Voy 



ix^ - 2x)(cot ^V^r/ - W^S) x^ + -X R^^ 

\ 6 3 3 / 



5x - -x^ - x^\ V^f], - ( Y - a; J Rrjt^ 



(75) 



with the BC 



ui = {x = 0). 



(76) 



It is easy to verify that 



/ X 



X 



wv\ - cot ev^T] 



X 



X 



X X 



-(To-^^r r""^[-6^T-rr"' 



X 

12 



X 



R^ 



•tz 



(77) 



is the solution. We note that the complete set of VCs (|A1CI|) obtained in Appendix ^ guar- 
antees that all the terms discarded in obtaining solvable ODEs for wi and Ui are small in 
comparison with (the biggest of) those terms that are retained. 

At this point, we could proceed to solve the x-NS and y-NS equations for the pressure 
and velocity corrections pi and vi respectively. However, these corrections are not needed for 
obtaining the second-iteration EE. (We will only need the pressure and velocity corrections 
for later iteration stages, and we calculate those corrections in Appendix p.) 
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IV. THE DISPERSIVE-DISSIPATIVE EVOLUTION EQUATION 



The (exact) kinematic condition (^) at x = h can be written in the form 



Vt + {vo + vi)r]y + {wn + wq + wi + W2 - 2)r]^ 



Uq + Ui + U2 {X = l + T]) 



(78) 



where we have used V = 2 (pB]). Dropping the terms containing unknown velocities (those 
with tilde) and using the Taylor series to relate the velocity components at x = h to those 
at a; = 1, we have 



Vt + {u;NxV + U!o)r]z = Uo:,;T] + Ui (x=l). 



(79) 



In (^), we have dropped the terms Vorjy and Wirjz as they are smaller than ui [see Eq. ([r5|)]. 
Also, recall that wjyxix = 1) = [see Eq. ([T^)]. Using the expressions (|36D, (0), and (|77[) 
for Wo, Uq, and Mi, we obtain 



V - j^RVz 



2 2 

+ "^VVz + -5ri-,-, - - cot er]yy 
t o 6 



where by definition 



+ -WV^T] + 2V'^ri, = 



(80) 



6 = -R-cote. 

5 



(81) 



However, the term oc Rrjtz ~ {R/Z)rit <^ rjt since R/Z <^ 1 [see ( [Al(j| )]. Dropping this small 
term, we have 



7]t + A7]7]z + -5r]zz - - cot 9r]yy 



+ -WV'r] + 2V'r/, = 0. 



(82) 
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Simple linear- stability analysis can reveal the dynamical role of some terms here. As- 
suming an infinitesimally small disturbance in the form of a normal mode, r] oc exp(st — 
iwt) expi{jy + kz), it readily follows from the linearized version of ( |82D that 

s = - [6^ - (cot e)f - W{k^ + 2fe + f)] (83) 

o 

and 

uj = -2k{k'^+f). (84) 

Here s is the growth (or decay) rate for the disturbance and so the third, fourth, and 
fifth terms in Eq. (|82|) which give rise to growth (or decay) are dissipative [considering the 
destabilizing term (the one with 6) as a negative dissipation]. In contrast, the last term in 
(^) only makes a contribution to the (real) frequency u, rather than to the growth rate s, 
i.e. it does not lead to growth or decay of disturbances. Thus, this (third-derivative) term is 
dispersive. 

Clearly, for instability to develop (i.e. for s > 0), we need 

(5 > 0, (85) 

a condition we assume fulfilled from now on. This yields the so-called critical value Re of 
the Reynolds number, R^ = (4/5) cot 9, at which the instability sets in. 

One can see from the above derivation of the dispersive-dissipative EE (^) that the desta- 
bilizing (third) term originates from the inertia terms of the NS equations. The (stabilizing) 
fourth and fifth terms are due respectively to hydrostatic and capillary (i.e. surface-tension) 
parts of the pressure. Finally, the last, odd-derivative term is due to the viscous part of the 
pressure. Such a purely dispersive term also appeared in the EE obtained by Topper and 
Kawaharali3 for an almost vertical plane: they used the small angle of the plane with the ver- 
tical as their (single) perturbation parameter (see also the discussion in the Introduction of 
the present paper). Our derivation shows that assumption to be unnecessary. In particular, 
for the vertical film cot ^ = 0, and Eq. ( p2|) becomes 
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Vt + 4^VVz + Y^RVzz + -WV'r] + 2V\ = 0. (86) 

Although an equation of this structure (but with arbitrary coefficients) was postulated as 
a model equation in Ref. p5|, it cannot be obtained from the derivation of Topper et al.": 



Since their small parameter is proportional to cot 6, it becomes zero for the vertical case, 
and the Reynolds number (also proportional to the small parameter in that SP derivation) 
vanishes — which, clearly, cannot correspond to any flow at all. 

The (infinite-dimensional) dynamical system governed by the dissipative equation (|82D 
essentially forgets initial conditions as it evolves towards an attractor. There may be fluctua- 
tions on the attractor, but there is no systematic change in time. Clearly then the amplitude- 
decreasing, stabilizing term must balance the destabilizing one (the latter tends to increase 
the deviation amplitude). So the two dissipative terms are necessarily of the same order of 
magnitude on the attractor. 

As to the magnitude of the dispersive term relative to that of the dissipative terms, 
there can occur, depending on location in the parameter space, each of the following three 
possibilities: (1) these terms are of the same Om] (ii) the dispersive term is small (and then 
the amplitude is determined by the balance of the nonlinear term with the dissipative terms); 
and (iii) the dissipative terms are small (and the nonlinear term balances the dispersive one). 
In the ffist case, continuing the iteration process would lead to small corrections to the terms 
which cannot significantly change the evolution. In the second case, the small dispersive term 
can be omitted with a negligible effect, so the corrections would be again immaterial. 

But in the third case, when dissipation is small, the situation is very different. Discarding 
the dissipative terms leads to a 2D Korteweg-deVries (KdV) equation which was simulated 



numerically in Ref. ^ The KdV equation is purely dispersive and never forgets the initial 
conditions. It has a one-parameter family of axisymmetric solutions which are traveling 
solitons, similar to the well-known ID KdV case. Depending on the initial state, there may 
be solitons of different lengthscales (and therefore moving with different speeds) in the final 
state. However, if the small dissipative terms are present, they will slowly change the initial 
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soliton of an arbitrary lengthscale. It will evolve along the soliton family until the lengthscale 
is attained which provides for the balance between the two dissipative terms (this effect was 
first studied for the ID case in Ref. ^; see also Refs. |28|^0D. Therefore, the dissipative 
terms, even when small, are important: they determine the lengthscale of the solution. 

However, only the largest-magnitude terms are guaranteed to be correct in the above 
beginning-iteration derivation; as for the smaller terms, further iterations might yield sig- 
nificant corrections to them. We consider this question (of higher iterations and corrections 
to small dissipative terms) in the next section. It turns out (perhaps, surprisingly) that 
such corrections can be important locally under some parametric conditions, but that no 



(single) corrected EE can approximate the evolution for all time. Equation (|82|) is thus the 
most general of those EEs that can be valid globally in time — under appropriate parametric 
restrictions, which can be completely determined only with the analysis of higher iterations 
of the NS problem, as is done in the next section. For the rest of this section, we continue 
the consideration of the GEE (^2]). 

From the condition 6 > 0, Eq. (BH), it follows that 



5 

i? > - cot ^ > cot e. (87) 

4 

From Eqs. (J8^) and (|83|) , it is clear that in order for instability to develop, we need Srjzz > 

{cot 6)riyy. Using the y and z lengthscales, this yields, Y^/Z"^ > {cot9)/6. Noting that either 

(5 -C -R ~ -Re ~ cot 6 or 6 ^ R, we see that Y ^ Z or Y r^ Z, except perhaps for R ^ cot 6. 

For simplicity, we assume that Z < Y, which seems to be the case in all experiments we 

know about. Then 

L = min(Z, Y) = Z. (88) 

With this and the condition (p7[), the VCs ( |A15| ) reduce to 

1 R W 



max 



'^' L2' L' L3 



< 1. (89) 



These are conditions of instantaneous validity; they involve the local (i.e. instantaneous) 
parameters A{t) and L{t). 
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As was discussed above, due to the dissipativeness of the EE (0), the system evolves 
towards an attractor, and in the asymptotic hmit of large times we have A{t) = const = Aa 
and L{t) = const = La. Since (similar to Refs. 0,|l3) the destabilizing inertia term should 
be of the same Om as the stabilizing, capillary one, i.e. Srjzzi^ ^^/-^a) ~ WV^rj^r^ WA/Lf), 
the (dimensionless) characteristic lengthscale at large times La can be taken to be 

i.= (y) . (90) 

Similarly, the asymptotic magnitude of the characteristic amplitude Aa is determined by 
the balance between the nonlinear "advective" term and either the dispersive term or the 
capillary one (whichever is larger): Aa = m.a.x{W/Ll, l/Ll). Using these asymptotic values 
of parameters, the conditions (|89D can be written as m.ax(W/Ll^, R/La, l/-^a) <^ 1- Noting 
that in view of Eq. (|0D, W/Ll = 5 /La and [see Eq. (|l|)] R = (5/4) (5 + cot 9) > S, so 
that W/L^^ < R/La, we can simplify the VC to max {R/La, L^"^) ^ 1; in terms of the basic 
parameters, 

„.J,^A«i. ,.£^fl(A)"%<i. (91) 

In the next section, it is shown that we also need 

maxfi?, R^) , . 

^" V «^ (""' 

(otherwise, the dissipative terms contributed by higher iterations can become significant, 
and the evolution cannot be all-time-describable by a single EE). All three parameters, a, 
(5, and 7, are small if (recall that 6 < R) 

R pq'^h^sm'^9 _ ,„„, 

and 

So, if we are in the domain of the space of basic parameters which satisfies the condition 
max(aR,/3ij) ^ 1 [or a bit more general, but less simple condition max(a,/3,7) ^ 1], then 

27 



Eq. (|8^) is good for all time [provided the initial amplitude and lengthscale satisfy the 
conditions (|^)]. Therefore, we call such conditions the "global" VCs. 

We can transform the GEE (^) to a "canonical" form — which contains only two "tun- 
able" constants — by rescaling the variables with appropriate units: 

1] = Nfj, Z = LaZ, 

y = Lay, and t = Tat, (95) 

where A^ = 1/(2L^) and Ta = L\/2. Dropping the tildes in the notations of variables, the 
resulting canonical EE is 

m + VVz - i^Vyy + ^'^Vz + e [t]^^ + V^T]^ = 0. (96) 

The control parameters in this equation are 

1 



-VW6 (97) 

o 



and 

L„cot9 1 IW 



^^^^^.,^^^^^ ^gg^ 



For the case of k = (e.g. flow down a vertical wall), Eq. (|96D becomes essentially the 
model equation postulated and numerically studied in Ref. ^. If, in addition, e = k = 0, we 
have the 2D KdV equation^, whose ID (dy = 0) version is the usual KdV equation. When 
K = but e — > oo (so that, after an appropriate rescaling, the dispersive term disappears from 
the equation), Eq. (^) becomes the one obtained by Nepomnyashchya, whose ID version 
is the Kuramoto-Sivashinsky equationE^l^a. AH these equations are thus limiting cases of EE 

Numerical simulationsc3 of this equation have shown that its solutions remain bounded 
for all time (which property always is implicitly assumed in global-validity considerations, 
and therefore must be directly verified after the EE has been obtained). For example. Fig. |l| 
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shows evolution of the surface deviation "energy" / rfdydz from an initial small-amplitude 
[t] ~ 10^^) "white-noise" condition to a statistically steady state (for k = and e = 1/50). 
Detailed numerical studies of Eq. (^) will be presented in a sequel to this paper. 

V. ADDITIONAL DISSIPATIVE TERMS 

In this section, we examine the implications of the additional dissipative terms arising 
from further iterations. We take into account explicitly all the linear dissipative and disper- 
sive terms with derivatives of order four or less by going through the iterative process twice. 
We also estimate the effect of dissipative terms with derivatives of order six or more on the 
EE. 

We note that for obtaining the evolution equation, we need only the successive iterates 
of the normal velocity u. This is because the nonlinear terms involving wrjz and vrjy in the 
kinematic condition make smaller contributions to the EE. (Indeed u ^ w^ + Vy, and e.g. 
Wzi] generates the same terms as the Wz part of u, but with the extra small factor rj.) We 
have found the next two iterative corrections for the normal velocity, U2 and M3, Eqs. ( [B22|) 
and (|B37| ) in Appendix p. [We need U3 (in addition to U2) because its dissipative terms are 



not guaranteed to be much smaller than those of U2- However, the dissipative terms of u^ 
are much smaller than those of U2, so we do not have to consider U4.] Using these in the 
kinematic condition (^), we obtain [see Appendix ^ for details] 



Vt 



—Rriz R cot eV^ri + i?V 

12 ' 15 ' 672 ' 



2. 2 „ _o 2. 



^mz-^5Vzz+^ cot er]yy-2V^r]z--WV\ 
too 6 



23 

— i?-2cot0 

15 



{rj 



5 „ .^2 2 „2 6 ^. 331 „„2 1241483 . 477523 ^o n^2 

H Rcot eV^riz - -R Vzzz + - cot 9V^n RV^Vzz R Vzzzz H R cot 9V^nz 

14 '7 ' 5 ' 168 ' 8108100 ' 2494800 ' 

(99) 

We get rid of time-derivatives in RHS by twice iterating this equation, substituting (for each 
time-derivative in the RHS) the RHS (^) of (|99D itself; the remaining time-derivatives in 
the final RHS are omitted, because further iterations would only lead to derivatives of an 

29 



order higher than four, with small resulting contributions. [In fact, in the first iteration, 
when substituting into the (mixed-derivative) terms with the second spatial derivatives, it 
is enough to retain from the RHS (p9| ) only the terms with no time-derivatives, and with 
space-derivatives of the order two only. As for the term with the first spatial derivative, 
oc ritz, the terms without time-derivatives and with space-derivatives of orders two and three 
are substituted into it, and also the term {5/12)Rritz itself, with the result {5 / 12)"^ R'^rjizz- 
In the second iteration, it is sufficient to only substitute into the latter term, and — since 
it already contains two spatial differentiations — only the purely spatial second-derivative 
terms should be substituted into it.] As a result, we obtain the following equation: 

2 2 40 40 2 2 157 

Vt+4:VVz+^5rizz-- cot eriyy+2V^riz+—RSrizzz—^Rcoterizyy+-WV^r]-- cot 9\/^ri+—R\/^riz 

8^ 2.^4 1213952 „ 138904^2 ^ „ /16„ ^ n\ ^ ^ 

+ -R cot^ ev^, + ^^^^^R%... - j^^^R' cot ev%z +{-R-2cot e) (,,.). = o. 

(100) 

The ID {dy = 0) limit of this equation coincides with the small-amplitude limit of the 
EE obtained, with the same numerical coefficients, in Ref. 0, but our 2D version is new. 
(We remark that there are several mistakes in the presentation of steps leading to the final 
equation, Eq. (27) in Ref. 0. However, Eq. (27) itself appears to be correct. The same 
numerical coefficients appeared in an even earlier paperHil in a linearized ID context.) The 



(2D) terms with derivatives of order three or less agree with Ref. |T5|, and those plus the 



surface-tension (W) term — with Ref. 0. However, note that some terms of (|100|) are always 



negligible. For example, the two dispersive third-derivative terms containing R are clearly 
smaller than the corresponding second-order dissipative terms (because of the instantaneous 
VCs R/L ^ 1 and cot6/L ^ 1), and therefore are negligible in all cases. 

As was mentioned above, since we are only interested in situations where persistent 
nonlinear waves are present, 5 > 0, we have either 6 ^ R or 6 <^ R (but not 6 ^ R; 
see Eq. (^Tj). When 6 ^ R, the additional fourth-derivative terms in ( P^ ) are each smaller 
than the destabilizing second- derivative term (because of R/L ^ 1 and/or 1/L^ ^ 1). If 
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max(i?, R^) <^ W [see Eq. (p2[)], those terms are much smaller than the stabilizing capillary 
term, and we return to the GEE ( p2D with global VCs (p3| ) and ( p^ (this holds as well even for 
6 <ti R). But if W is not large enough, so that Eq. (|9^ ) is violated and the capillary fourth- 
derivative term is much smaller than (at least) one of the non-capillary fourth-derivative 
members, then the destabilizing term cannot be balanced, and the EE leads to the unlimited 
growth of amplitude. Thus, clearly, at such parameter-space locations the EE ( p. OOP can be 
valid locally, i.e. for a limited time only, but it is not globally-good. 

Consider now (for the rest of this section) the complementary case 5 -C i? (so R ~ Re), 
and with W sufficiently small, so that the condition max(i?, R^) <^W is violated. We have 
obtained, except for numerical coefficients, all the essential terms (which turn out to be 
linear; see the end of Appendix jC]). Including these terms in ( p.00| ), the EE can be written 
in the general form 

2 8 2 

Vt + ^VVz + -Sr]zz - -TrRcVyy + "^Vzzz + -Wr]^^^^ 



+ -Rc{vVz)2 
5 



2581 32 ■ 



1400"^' 3378375 



+ sum of terms of type Rf+^r]'^^^^ = 0. (101) 

Here the superscript on rj, enclosed in parentheses, refers to the order of the spatial derivative; 
/ > A; > 0; A; = 0, 1, 2, ■ ■ ■; / = (A; + 1), (A; + 2), ■ ■ ■; and every {k + /) is odd. In ([TOTD , we have 
used the leading Taylor-series approximation putting R = R^ (recall also cot6' ^ 4_Rc/5). 
Also, F ^ Z (as a consequence of 6rizz > RcHyyi which is required for instability to develop). 
Furthermore, we have not included additional nonlinear or dispersive terms in Eq. ( |101| ): the 
dissipative nonlinear terms can be shown (see Appendix |C|) to be smaller than the leading 
nonlinear dissipative term oc R{7]7]z)z, and all the dispersive terms are smaller than the linear 
one V'^rjz. 

Normally, the coefficients of the terms are Om(1)- However, sometimes a coefficient is <^ 1 
because of an accidental near cancellation of terms, as e.g. is the case for the term oc R^rj^zzz 
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in ( |101|) . Then we say the term is "degenerate" . Comparing the additional dissipative terms 
with the term oc R^rizzzz and taking into account the instantaneous VCs R/L <^ 1 and 
1/L^ <^ 1, the only terms which may not be negligible are those of the structure R^^~^ri^'^^'> 
(n = 2, 3, ■ ■ ■) (even those could have been neglected in comparison with the term oc R^rjz^zz 
were the latter non-degenerate). Hence, the EE can be simplified: 

2 8 

Vt + 4:VVz + -Srizz - T^RcVyy + '^Vzzz 



+ -Rc{vVz)z + 
5 



2581 32 ■ 



1400"^' 3378375 



+ 5^c„i?f~V2") = 0, (102) 

Tl=3 

where c„ are numerical coefficients. Can this equation be valid for all time? For this to be 
the case, the destabilizing term (2/3)5?7^2 has to be balanced by a stabilizing one — by the 
term oc Rcf]zzzz (the other, degenerate fourth-derivative term is clearly destabilizing) or by 
the first non-degenerate and stabilizing higher- derivative term of (|102D , whichever term is 
dominant. 

Suppose first that the fourth-derivative term is the dominant stabilizing one. Then Eq. 



11021) becomes 



2 8 

r]t + 4:r]riz + -drj^z - -TzRcflyy + '^Vzzz 

8 2^81 

+ -Rcivvz)z + Y^RcVzzzz = 0. (103) 

The balance of the destabilizing term Srj^z with the stabilizing term Rrj^zzz yields the 
lengthscale L ~ Rc/S. Comparing the dissipative term with the dispersive term, we have 
RcVzzzz/Vzzz ~ Rc/L <^ 1: the dispersive term is always dominant. Therefore, it is the dis- 
persive term which is to balance the nonlinear one, which yields the characteristic amplitude 
~ 1/L^. Using this, it is easy to see that the nonlinear dissipative term ~ Rc{vVz)z is exactly 
Om of the linear dissipative terms; therefore, the nonlinear term plays a significant role. It 
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is destabilizing, and numerical simulations indicate that the solutions of ( |103|) blow up [see 
Fig. ^, which shows the blow-up of energy (in contrast to Fig. |1|) governed by the equation 
rescaled to the canonical form similar to Eq. (p6|), 



Vt + VVz + Vzzz - l^Vyy 

+ e [vzz + Vzzzz + (1120/2581) (r/r^,),] = 0, (104) 

with K = and e = 1; we have observed similar blow-up behavior with all values of e we 
tested in the range 10~^-10^]. The thickness of the real film, of course, is bounded uniformly 
for all time. Hence, Eq. ( |103|) (which can be good for a limited time) is not globally-valid. 
Physically, we believe that the growth of amplitude will be arrested by viscosity after small 
lengthscales develop, which would violate the sraall-R/L VC for the single-EE description. 
So the EE ( |103| ) cannot be valid for large times. 

It remains to consider the hypothetical case when a higher-derivative term oc R^'^r]^'^'^^ 
is the dominant stabilizing one. [In particular, this implies that ?> < n = m (say), c„ <^ 1 
for n < m, Cm ^ ^, and c^ has an appropriate sign: c^ > if ?7i is even and c^ < if 
m is odd.] We do not believe this actually happens: such terms are traced back to be due 
to the inertia terms in the momentum NS equations, the same inertia terms which give 
the destabilizing second-derivative term of Eq. (|102|) , and it seems unlikely that the same 
physical cause can be responsible for both stabilization and destabilization. Therefore, we 
believe the first nondegenerate term will turn out to be destabilizing. 

Based on the (linear) studiesE3 of the Orr-Sommerfeld equation, we have computed the 
next coefficient, which shows that, albeit stabilizing, the sixth-derivative term is again de- 
generate: C3 = -16173184/1718663821875 ^ -0.941 x 10"^ (This is also remarkably close 
to the fourth-derivative coefficient, C2 ~ —0.947 x 10~^. If further c„ were all degenerate 
too and of the same order of magnitude, the destabilizing fourth-derivative term would be 
dominant, and hence stabilization and all-time valid EE impossible.). 

A special case to consider is when the destabilizing fourth-derivative term (which is 
clearly much greater than the sixth-derivative one) is nearly cancelled by the stabilizing 
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4th-derivative term aRj]^^^ (where a = 2581/1400). As the term csR^r]^^^ is stabihzing, 
oneE2 can ask whether there can be an all-time valid EE if Ic^Rlr]^^^ ^ \aRc — C2-Rj?)?7*^^'*|. 
Our answer to this question is as follows. In this case, aRci]^^^ ~ |c2|-Rc'7''^^ ^ \c3\Rlr]^^\ Since 
the dispersive (third-derivative) term must dominate the fourth-derivative term aRcrj^^\ it 
dominates even stronger the sixth-derivative term csRlr]^^^ [see the discussion following Eq. 



( |103| )]. As usual, the lengthscale is determined by a balance between the dominant linear 
dissipative terms, 5r]zz ~ Icsl-Rc^*-^-* <^ aRcXj^^^ ~ Rc{vVz)z [see the arguments following Eq. 



( |103| )]. Thus, the destabilizing nonlinear term is the greatest one; there is no other term 
which could serve as a counterbalance. We arrive at the conclusion that the hypothetical 
EE with the sixth-derivative term cannot be valid for all time. 

We have not attempted to determine the next coefficient, C4, because of the large volume 
of calculations which would be required. Based on what we have said above, we expect the 
eighth-derivative term to be nondegenerate and destabilizing, and then no single EE can be 



globally valid under the circumstances. It follows that the GEE (|8^ is the most general one. 
The validity condition W ^ R^ (which would be sufficient even if C2 were nondegenerate) can 
be relaxed a bit: it is enough to require that the capillary term dominates the (presumably 
nondegenerate) eighth-derivative one, W/L^ > Rl/L^- (With L^ ~ W/5 [see Eq. (pO])], we 
get W'^/'^/{5^/'^Rl) > 1.) The EE ( p^ ) is valid locally only, under the instantaneous VCs 

If, however, the first nondegenerate term oc i?^"-"^?]*^^") turned out to be stabilizing, the 
EE ( |102D would have a domain of global validity, albeit a very limited one. It is straightfor- 
ward to find the corresponding global VCs. Indeed, the balance between this term and the 
destabilizing one, /?2m-i^(2m) ^ Srfzzi yields the lengthscale L, 

D2m— 1 



-2m-2 



(105) 



Using this lengthscale, the "modified- i?" VC takes the form 

< 1, (106) 



/ 5 \ ^/(^™-^) 



L \RcJ 
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and the small-slope VC becomes 

1 

L2 



l/(m-l) 

< 1. (107) 



The (relaxed) conditions of dominance are Wrj^'^'^ <^ R^^'^rj^'^'^^ and RcTj^^^ <^ Rf^~^ri'^'^"^\ 
i.e. W < i?2m-Y^2m-4 ^^^^ R^Rjlf"^-^ > 1 (which with ([T05|) are easy to recast in terms 
of the global parameters only). The ratio of the nonlinear dissipative term, oc Rc{vVz)z, to the 
stabilizing one, oc R^J^~^ri^'^'^\ is A/{Rc/ LY™'~'^^ which is required to be small — otherwise, 
we can have the blow-up of solutions, similar to the case with the term oc Rci]^'^'' being 
dominant. The comparison of dissipative term, oc Rl"^~^ri^'^"^\ to the dispersive term, oc rj^zz, 
shows that the latter can be greater or smaller than the former depending upon whether 

When i?^™"^///^™'"^ <^ 1 (dispersion is large), the characteristic amplitude, obtained by 
balancing the term rjrjz with the dispersive term, is l/L^ [hence, the small-amplitude VC, 
A <^ 1, coincides with ( flOT] )]. Otherwise, i.e. if R^J^~^ / L'^'"^''^ > 1, the amplitude is deter- 
mined by balancing the term oc rjrjz with the dominant stabilizing term, A ~ (i?c/-^)^™~^- 
Using the global-parameter estimate of the lengthscale, Eq. ( |105| ), all the above conditions 
are readily reduced to certain global VCs (expressed in terms of global parameters only); we 
do not write them here, in view of the likely nonreality of the imaginary case of dominance 
of the term oc R^J^~^rj^'^'^\ 

VI. SUMMARY 

We have considered flow of a liquid film down an inclined plane. We have posed and 
studied the following questions: What are the least restrictive parametric conditions for 
which the wavy film flow can be approximated for all time by a single (local) evolution 
equation, and what (if any) is the most general form of such an equation? 

We have argued that the dissipative-dispersive evolution equation (p2D [which we derived 
by an iterative perturbation method of a multiparameter type] is such a general EE. Any 
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all-time valid EE derived by a single-parameter technique is necessarily nothing else but 
essentially the general EE in which some terms have been omitted. Also, the domain of 
validity of such a "partial" EE is necessarily a suhdomain of the "umbrella" domain of 
global validity given by Eqs. (p3| ) and (0). [In particular, in such domains the amplitude of 
waves is necessarily much smaller than the mean film thickness.] 

It is clear that any evolution equation which follows from a multiparametric approach 
(such as the iterative technique we have employed in this paper) can be also obtained with 
the conventional SP approach. However, the significant advantage of the MP derivation is 
that it covers at once all possible SP derivations of EEs (the number of which is, in principle, 
infinite in the SP approach, corresponding to the different choices of the small-parameter 
powers for the system parameters). Also, comparing the two derivations of even a particular 
EE, the MP derivation is justifiable for much less restrictive domains of the parameter space. 

The theory also yields the explicit approximate expressions in terms of film thickness 
for the pressure and components of velocity [Eqs. ([50|), (^), (^), and (^)], and thus a 
complete description of film dynamics. 

We have derived an EE ( |100D containing additional high-derivative terms which can be 
essential near the threshold of instability, R ~ Re-, if the surface tension is sufficiently small. 
However, this EE ( |102|) is good for a limited time only and cannot be good for all time. [The 



conditions of such a local (in time) validity are given by Eq. (|89|) .] Under certain parametric 
conditions, the (numerical) solutions of that EE blow up due to a nonlinear (quadratic, 
second-derivative) destabilizing term. 

The EE (^2]) is relatively easy for numerical simulations of the 3D waves in the inclined 
film. We have obtained good agreement with transient states and transitions observed in the 
physical experiments of Ref. ^. Under certain parametric conditions for which the dissipative 
terms of the EE are small, we observed self-organization (from the initial white-noise small- 
amplitude conditions) of unusual highly-ordered patterns of soliton-like structures on the film 
surface (the pattern consists of two travelling-wave subpatterns which move with different 
velocities). The studies of the evolution equation (|^) will be published elsewhere (Ref. |39|; 
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see also Refs. |,0). 

A similar analysis leads to analogous single-EE theory of a film flowing down a vertical 
"fiber" . We believe that similar theories can be useful for a variety of other systems. 
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APPENDIX A: VALIDITY CONDITIONS 

From (^), the estimate of Wi, in terms of the basic parameters and length- and 
amplitude-scales, is 

/I 1 R cote 

Wi ~ max — -, — -, — , , 

R W W 



T' Z3' zy2 

Estimating ui from (^), we have 

/I 1 i? cot 6^ cot 6* 
Ml ~ max 



A (Al) 



R W W W 



TZ' Z4' ^2^2' y4 



A. (A2) 



In obtaining a solvable equation for the boundary condition on wix at x = 1 [Eq. (|7TD], we 
have dropped the term Uiz{x = 1). This implies 

uiz < wix (x = 1). (A3) 

Using the Om estimates for ui, Eq. ( [A2D , and Wix{x = 1), Eq. (|7lD , the above requirement 
reduces to 
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(I 1 i? cote cote i? 

u^z\X = 1) ~ max — -, -, — -, — --, — --, , 

i-z\ ) \ Z^ ZY^ Z^ Z^ Y^ TZ 



W W W\fA\ . ^, A 

This again yields the conditions (^), p9|), (pOf), (|3lD, and, in addition, the following condi- 
tions: 

cote _, , _, 

^— < 1, (A5) 



Z\ fcote\ ,, , 

^ « 1, (A7) 



W 

< 1, (A8) 



ZY 
and 

(f)(^)-- <-' 

These conditions, along with (|TT]), (^Sj), (P^D, (PO]), and (0), form the complete set of validity 
conditions for the present theory. They are sufficient to justify all the simplifications of the 
equations. Thus, the complete set of VCs is 



max 



' 1 1 i? cot e z cot e i? 

A 

' ^2' y2 ' 2^' Z ' Y'^ ' T' 
W W (Z\ WZ 



^3' y2^' yy J Y^ 

Using these conditions, it is easy to see that wi <^ wq and ui <^ uq. 
From ( p6|) , we can estimate the Om (Rpi)'- 

A 1 1 i? cote cote 



< 1. (AlO) 



Rpi ~ max — , — - 



Z' Z3' ZY^' Z2' Z2 ' F2 ' 
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R W_ W W_' 
XZ^ Z^' 'z^Y^^ Y^ 



(All) 



Using the conditions (|A10|) , it is easy to show that Rpi/Z <^ Rpo/Z or pi <^ po- From 
( P12|) , we estimate the Om {vi) as 

A 1 1 R cote 



Vi ~ max 



YZ' YZ^' ZY^' YZ^' YZ^' 



cote R W W W A cote 
Y^ ' TY Z^ Y Z^^ Y^Z"^^ Y^^ Y ' 

AW AW Rcote RW 



YZ"^' ys ' ZY ' YZ^' 



RW Rcote 
ZF3' TY 



RW RW 
' TY'^Z^ TY^ 



A. 



(A12) 



Using the estimate of Vq [Eq. (^)] and the conditions ( |A10| ), it is easy to see that Vi -^ Vq. 
The VCs ( lAlOP guarantee that all the terms involving wq, uq, ^"0, po, wi, ui, vi, and pi that 
were dropped in obtaining solvable ODEs for the same quantities are small in comparison 
with the terms that were retained. 

Estimating the Om of various terms in Eq. (|8^) , and noting that 6 <^ R or 6 ^ R, we 
find that 



— ~ max 
T 



1 1 5 cot 6* Z cot 6* 
' Z2' F2' z' ~Z~' F2 ' 



w w wz 



^3' ZY^' Y^ 



and, consequently, taking (|A10| ) into account 



R 

— ~ max 
T 



1 1 (5 cote cotez 
' z2' F2' z' ^~' r2 ' 



(A13) 



W WZ w 



Z3' F4 ' ^y2 






(A14) 
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Hence, the parameter R/T is small as a consequence of the smallness of other parameters 
in ( [Aiq ), and thus can be omitted from there. Also, ( |87D yields cot 6 /Z < R/Z ^ 1, so that 



the parameter cot 6/Z can be omitted as well. The somewhat simplified validity conditions 
are 

1 1 R cotez 

max 



A 



w w wz 



< 1 (A15) 



Z3' V^Z' Y^ 
(which, in turn, significantly simplify [see Eq. (^)] iiY > Z). 

APPENDIX B: SECOND AND THIRD ITERATIVE CORRECTIONS 

To find the second and third corrections for the velocity components and pressure, we 
need first to determine pi and ^i, the approximations to the exact solutions pi and vi. 
Considering the x— NS equation, 

1 1 

Plx = -POx + ^(Mo + Ui):rx + ^V^(Mo + Ml) 

-{Uo + Ui)t - {uq + Ui){uq + Ui)^ - Vo{Uo + Ui)y 

-{wn + Wo + Wi){Uo + Ui)^ + 2(mo + Ml)^ 

+ [terms containing Vi, W2, or U2] (Bl) 

where 

U2 = Ui — Ui. (B2) 

Using our standard procedure, the simplified equation is 

1 1 . 

Plx = -f^Uixx + 75V Mo - WnUqz + 2uoz - Moi 
K K 
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-(x - l){WV^r] - cot eV^T]) + r]t, 
ri 



x^ - -x^ + -) Vzz + (a;' + 2x- 5)-^ 



(B3) 



where Eqs. (^) and (|3I1| ) for pjv and po have been utihzed. The normal-stress balance con- 
dition 3X X = h yields 
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'[{uo + Ml), + t/7ijr7j {l+ryj+r^2|-i_^ J^^^ |^^ ^ ^2^ ^ ^^^ |^^ ^ ^2^ _ 2r^^r/,r/,J (l+ryJ+r^2)-3/2 



-|- [terms containing Wi, W2, or M2] (a^ = ^) 



(B4) 



which simplifies to the BC 



Pi = -POxV + —UoxxV + -TzUia 

K K 



ri o 



20^2 2 4 , ,, 

3R ^' ~ R^^' ^ 3^*^ (a; = 1). 



(B5) 



The solution is 



Rp^ = (a;2 _ 2x - 1) WV^rj - cot 0V^r/ 



fx^ x^ 4 7 \ „ 



a; 



+ ( — + a;2-5x-3j V^r^, 



+ {x + -] Rr]tz - 2r]r]z 



(B6) 
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The I/— NS equation is 

Vlxx = -Voxx + R{PN +Po+ Pl)y - '^'^Vq + RVqi 

+R{uq + Ui)vo^ + RvoVoy + R{wn + Wq + Wi - 2)t;o^ 

+ [terms containing W2, U2, P2, or Vi], (B7) 

where pi = P1+P2, and vq is given by (|58D . The simphfied equation (for the approximation 

vi) is 

Vlxx = Rpiy - V^t;o + R{wn - 2)t;o^ + Rvot 
or, in terms of 77, 



Vlxx = -(2x2 - 4x - 1) cot eV^riy + —x^ + -x^ - —x'^ + lOx^ - v^ " T7 ^^^'v 



2 5 2 4 16 3 ^^ 2 28 7 

— x^ + -X x^ + lOx X 

15 3 3 3 15 



I X 1\ 

-{x^ - 4x^ + Gx^ - 4:x)Rcot9r]y^ + (x^ - 2x)Rcot9rjyt h x^ - 5x Rrj^yt 



+ Pi{mz)y + P2W^V^r/, + P^RWV\, + P4i?W^V277,, + P5i?W^V'r?,i, (B8) 

where vi is the approximation to the exact solution vi. The coefficients of terms that are 
hkely to contribute to spatial derivatives of order five or more, as well as the coefficients 
of nonlinear terms which are small or dispersive in nature, are represented by Pi, • • ■ , P5, 
which are polynomials in x; we will not need their explicit form, just the fact that they 
are ~ 1. Henceforth, coefficients of type Pm will always refer to such polynomials in x. The 
boundary condition aX x = h is: 

vix = -Vox - {uq + Ui)y + {2[vQy - {uq + ui)^]r]y 

+ [(mo + Mi)2 + wix]r]yr]z 
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2l-l 



+ [voz + {wo + wi)y]r]^} [1 - r]y] 

+ [terms containing Vi, W2, or U2] {x = h), (B9) 

where we have taken into account (|37|). The simphfied BCs are 

Vlx = -VOxxV - UOyxV - Uiy - 2Uox?7y + Woyr]z 



2 5 8 

■-WW\ - —RVtzy + Y^RVzzy 



+Ci{r]T]z)y + C2(cot e)ririy + C3WriV\ 

+ C4WV\ + C,V\, (x = 1) (BIO) 

and 

^;^ = (x = 0) (Bll) 

where Ci, ■ ■ ■ , C5 are constants. (Henceforth, coefficients of type C^ are constants.) The 
solution is 

/x^ 2 o x2 ^ \ ^^2 / x^ x^ 4x5 5x^ 14 , Tx^ 158x\ „ 

vi = -{ x^ h 3x cot 9V^ri,,+ \ \ x^ \ Rn^zv 

I 6 3 2 / 'M 315 45 15 6 9 30 45 / ' ^ 



' x^ x^ x^ 2 o 4 \ „ , /x^ x^ 2x\ „ 

, 30 5 2 3 5 / '^ V12 3 3 / '^ 



1^ + ^ - ^a;=^ - y + 2x j i?r/t,^ + P6(w.)y + ^7(cot e)7]7]y + PsWr]V\ 



+ P,WV\ + PioV\z + PnRWV\, + P,2RWV\f (B12) 

By using the same procedure as above, we proceed to calculate u)2, ^2, P2, and 'i;2. The 
equation for W2 is 
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W2xx = -{wo + Wi)^^ + R{po + Pl)z - V^(wo + Wi) + R{wo + Wi)t 
+R{uo + Ui){wn + Wo + Wi)x + R{vo + Vi){wo + Wi)y 

+R{wn + Wo + wi - 2){wo + ^1)2 
+ [terms containing W2, U2, V2, or P2]. (B13) 



This yields 



W2xx = Rpiz - V^Wi + RWit + RUiWnx + R{wn - 2)^1^ 
+RUqWox + RWqWqz- 



Hence 



W2xx = (1 + 4a: - 2x^) cot dV% - f |- + |- - 1^ + 5a;2 + ^ j /^v't]. 



+ ( — - ^- + 12a;2 _ _ _ _ j Rr^^^^ - (a;^ _ Ax"" + Qx" - Ax) cot % 



x^ ?,x^ Ax^ Ax^ Ax^ 8a; \ n ( x^ 2x^ x^ 

10 5 3 3 3 3 / ' I 3 3 3 



+ [^ - ^ + 2a;2 j RcoteV^ri-l^ _^ + l^ + a;2_l^j R\^^jr^ -2x) R cot 9ritz+2x^RvVz 

+Pl3iVVz)z+Pl4WV%+Pi5V%z+Pl6Vzzzz+Pl7Vyyzz+Pl8RWV%z+Pl9RWV^ri+P2oRV% 
+ P2lRWV\z + P22R\t + P23RVzzt + P2ARVyyt- (B14) 

The boundary condition ai x = h is 

W2x = -{U!0 + U!i)^ - (Mo + Ul)z + {2[iwo + Wi)z - {Uq + Ui)^]r]z 

AA 



+ [(mo + ui)y + (fo + vi)^]r]yr]z 

+ [{vo + vi)z + {wo + wi)y]riy}{l - rilY^ 

+ [terms containing W2, U2, V2, or P2] [x = h). (B15) 

Dropping all the unknown and smaller terms, we have, at a; = 1, 

W2x = -WixxV - Uq^^T] - Uu + 2(Wo^ - Uox)riz + WQyT^y 
c o 

= —RVzzz - - cot eW^r], + 2Rriri, - 2 cot O'q'q, + C^Wr^y^'q, 
15 6 



+C7{riVz)z + C^{my)y + C^WV% + CioV'r/, 



2., 

zz 



+ CnRvtzz + CuRvtV (a; = l). (B16) 

Finally, the no-slip condition is 

W2 = (x = 0). (B17) 



Integrating (pi4|) with the boundary conditions ( pi6|) and ( |B17| ), we obtain 



^ V 6 3 2 J '^ V 126 60 60 12 9 180 y ' 

/x^ 17x^ 4 ITx^ 7x2 721x\ / x^ x^ 2x^ x^ x^ 4x=^ 232x\ . 

+ h X* ^ RVzzz+ \ \ \ \ R rizz 

V90 60 9 30 180 J ' \ 560 70 45 15 9 9 315 / ' 

/ x^ x^ x^ 2x^ 4x\ ^ „ / x'^ x^ x^ x^ x^ 89x\ „ 

V 30 5 2 3 5] ' V126 45 60 6 3 180/ '™ 

( x^ 2x^ x^ 2x\ ^„, / x'^ X® x^ x^ 5x3 199x\ , 
+ 77--^ + — -TT RcoteV^T] + -— - + 77-77-777 + ^^ TTTT i?V 



2 



^45 15 6 15/ 'V 252 45 15 12 9 180 / 

Cx^ x^ 2x\ /x^ 4x\ , , r, 

— - y + y 1 i? cot %, + I y + y 1 ^^^^ ~ 2x cot 0r/r/, + P25{vVz)z + P26W7]V\ 

45 



+P27{VVy)y+P28WV%+P29V\,+P3QVzzzz+P3lVyyzz+P32RWV\,+P33RWV^r]+PuRWV\, 
+ P35R\t + P36RVzzt + P37RVyyt + PssR'^^ + P39RVVt- (B18) 

The equation for U2, obtained from the incompressibihty condition, is 

U2x = -{Vl + V2)y - {W2 + W'i)^ (B19) 

where ^2 = ^2 + ^3. Dropping the unknown terms, we obtain 

U2x = -Viy - W2z 

or 

fx^ x^ x^ 2x^ Ux\ „ ,_^ / x^ x'^ 2x^ x^ x^ Ax^ 232a; \ „. 
^"^ \m 15 3 3 15 ; '" \hm 70 45 15 9 9 315/ " 

/x^ 2x3 a;2 \ ^ , ( x^ x^ x^ Ix^ V^x^ Ix^ 221a; \ ^ 

+ + 3x cot 0V^r7+ + + + R^ r]zz 

V6 3 2 ; ' V 315 60 30 12 9 30 45 / ' 

I x'' x^ 1x\ „ ^„9 / x'^ x^ x^ x^ hx^ 199a;\ ^9 /x^ 4x\ „, , 

\ i?cot6'VV+ \ \ \ Rritzz-\ \ R(vVz)z 

I 12 3 3 / ' 1 252 45 15 12 9 180 / ' I 6 3 / ^ ' '^ 



+2x cot 9{r]r]^) ^+P3g{rir]^) ZZ+P40W {r]V^r]^) ^+P4^i{r]r]y)y^+Pi2{vVz)yy+Pi3{VzVy)y+P4.4.WV'^r], 



4~, 
'zz 



+P45^^Vzzz + PmVzzzzz + PijTlyyzzz + Pi^RWV^T],,, + P^^RWV^ + PmWV\y + P5lV%yy 



-P52RWV%yy + P53RWV\yy + P^iRVtzyy + P55RWV\,, + P56R\tz + P^jRVtzzz 



+ P^f,R7]tyyz + P59R^\z + PmR{mt)z (B20) 

where U2 is the approximation to U2- The no-shp boundary condition is 

M2 = (x = 0). (B21) 

Integrating Eq. (|B20| ) with the no-shp EC, we obtain 
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X 



^ x^ 2x'^ x^ x^ x*^ llQx"^ 



V5040 560 315 90 45 9 315 / ' 



/ rfi I /yivJ rf'-J rf^ I 'y \ I 'TT''-^ rr*^ rf'~^ ~\'y \ 

+ + + LRcot eV^T], + + cot ev^ri 

1 630 90 15 6 15 / '" 1 30 6 6 2/ ' 



/ x^ x^ x^ Ix^ 19x^ 7x3 221x2 \ ^ 

+ \ \ \ R^ r\^z 

\ 2520 420 180 60 36 90 90 / ' 



/ x^ x^ x^ x^ 5x^ 199x2\ „, I x^ x^ x^\ ^ ^ /x^ 2x2\ ^, 

+ \ \ \ R Vtzz- \ LRcot Oritzz- \ RivVz 

1 2016 315 90 60 36 360 / ' 1 60 12 3/ ' 1 30 3 / ^" 



+x2cot6'(r7?72)^+P6l(W2)2^+-P62(Wy)?;^+-P63(W^)ra+-P64(^2^?y)?y+-P65W^('7V2r7^)2+P66M/^VV2 
+P67V\zz + PesVzzzzz + P69Vyyzzz + PjoRWV^,, + PnRWV\ + P72WV\y 



-P73V^Vzyy+P74RWV^7],yy+Pr5RWV\yy+P7eRVtzyy+P77RWV\zz+P78R\tz+P79RVtzzz 



+ PsoRVtyyz + P%lRy\z + P82R{VVt)z- (B22) 

In order to obtain all the fourth-order linear dissipative terms of the EE, we have to proceed 
to the next approximation [since (see Appendix 0) the dissipative terms of u^ are not 
guaranteed to be much smaller than those of U2', however, the dissipative terms of u^ are 
much smaller than those of M2]. Expressing U2 = U2 + u^, the incompressibility condition 
leads to an equation for the approximation u^: 

Usx = -V2y - W3^, (B23) 

(where V2 and w^ are the approximations to V2 and w^ respectively). The solvable ODEs 
for V2 and w^ will contain p2, the approximation to p2- Terms from p2 will not make any 
contribution to linear terms with derivatives of order four or less and the nonlinear terms in 
the final EE. However, in Appendix P, we will obtain an equation for the general iterate Pn, 
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and the form of the x— NS equation for p2 will guide us in formulating the general equation 
for pn. 

The x— NS equation for pressure p2 can be written in the form 

P2x = -{Po+Pl)x + -siUo+Ui+U2)xx + -s'^'^iUo+Ui+U2)-iUo+Ui+U2)t-iUo+Ui+U2)iUo+Ui+U2)a 
-{Vq + Vi){Uo + Ml + U2)y - {wn + Wo + Wi+W2~ 2){uo + Mi + U2) z 

+ [terms containing V2, W3, or u^ ] (B24) 



which yields 



U2xx V^-Ui 

P2x = -PixV + -p- + — Uit - 2Uu 



- WnUu - UqUqx - WqUoz- (B25) 

The boundary condition on pressure aX x = h is 

2 

P2 = -{pN+Po+Pl) + ^{{uo+Ui+U2)x + {vo+Vi)yr]y + {wo+Wi+W2)zVl-[{uo+Ui+U2)y + {vo+Vi)x]r]y 

+ [{vo + vi)z + {wo + wi+ W2)y]r]yr]z - [(mq + mi + ^2)2 + (wi + W2)x]Vz}{^ + vl + vl)~^ 
-a [r]yy (1 + Tit) + rizz (l + V^) -'^VyVzVyz] {I + V^ + vIT'"'' 

+ [terms containing V2, W3, or Us] {x = h). (B26) 



This yields, at a; = 1, 



2 



P2 = -PlxV + [^{UlxxV + U2x - {UOy + Vox)'ny 



-{uqz + wix)r]z] - (y[riyyrjl + rjzzill - 2r]yT]zT]yz] 
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-(3/2)aV'viVy + vi) (^ = 1). (B27) 

As mentioned before, p2 will not contribute to linear terms with derivatives of order 4 or 
less and the nonlinear terms. Hence, we will not solve for p2. 
Representing p2 as P2 = P2 + Ps, the y— NS equation for V2 is 

V2xx = -{vo+Vl)yy+R{po+Pl+P2)y-'^^{Vo+Vi)+R{Vo+Vi)t+Riuo+Ui+U2)ivo+Vi),^+R{vo+Vi){Vo+Vi)y 

+ R{wn + WQ + W1 + W2- 2)(t;o + Vi);^ + [terms containing V2, W3, U3, or p^] (B28) 
The simplified equation for the approximation f 2 is 

V2xx = Rp2y - V^Vi + Rvit 

+ RuqVq^ + R{wn - '2)vu + RwqVo^. (B29) 

The boundary condition on V2 a.t x = h is 

V2x = -{vo+Vi)^-{uo+Ui+U2)y + {2[{vo+Vi)y-{uo+Ui+U2)x]Vy+[{uo+Ui+U2)z + {wi+W2)x]VyVz 

+ [{vo + Vi)z + {wo + Wi + W2)y]riz}{l - Vy)^^ + [terms containing V2, W3, or u^] {x = h). 

(B30) 

The simplified EC on f 2 at x = 1 is 

V2 = -^VoxxxV'^ - VlxxV - {ijUQy^^Tf + Uiy^T] + U2y) 
+ [2{vq - UoxxV - Uici:)r]y + 

+ {voz + Woy^Tj + wiy)ri^] (x = 1). (B31) 

The no-slip condition on V2 is 

V2{x = 0) = 0. (B32) 
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As mentioned before, we will not solve the equation for V2 to find all the terms in it. We are 
only interested in the terms that make linear dissipative contributions of order four or less 
to the EE. Only the term involving R{wn — 2)f i^ makes such a contribution. There is no 
contribution from the tangential stress and the no-slip conditions. Thus, 

vi^x = R{wn - '2)viz 
~ \30 "Is'^^O 15~^~ 



Ix^ 8x^ ^x 



^^ — + — ] R' cot er]y,, (B33) 

where the superscript on V2xx indicates that RHS of the above expression contains only linear 
terms with spatial derivatives of order 3 or less with all other terms omitted. 
Solving for ^2 , we obtain 

(<3) _ / x^^ x^ 29a;S 31a;^ 7x^ 
^^ ~ V2700 ~ 270 ^ 1680 ~ "630" ^ "90" 

2x^ 2x^ 4x^ 344a; \ „o , ,^^,, 

+ R^ cot %^.. B34 

75 15 15 945 / ' ^ ^ ' 

Similarly, proceeding to W3, we find that only two terms in the z— NS equation make 
contributions to fourth derivative terms in ^3, i.e. 

W^x = RWNxU2 + R{wn - 2)w2z, 



or 



rf -'-^ <^y 



19a;S 47a;^ 2x^ Ax'' 8x* 8x^ 232a;2 464x\ 



(-<q') / X X LUX ^IX Z,X ^X OX OX ZOZX ^U^X \ „q 

3xx ^-^20 72 315 315 15 15 9 9 315 315 / ' 

^ %x^ QAx'^ 19a;6 IQx^ 2x^ 4x^ 4x^ 2x^ 4:x\ „2 n^2 



315 315 30 15 3 5 5 3 15 



+ ^- + — i?^ cot eV^T], 



x^ Ax'^ 29a;^ 31a;^ 7x^ sx" »x" ax 



3 fi^2 „ , . 

2 



+ + + + — i?^ cot dn,,,. (B35) 

V30 15 30 15 3 15 5 5 y ' ^ ' 

50 



We note that only the terms with R cot OV'^rjz and R^rjzzz in U2 [see Eq. ( |B22| )] make a 
contribution. Similarly, only the terms with RcotOrjzzi R'^Vzz, and RcotO'V'^r] in W2 [see Eq. 
( [B18| )] make the relevant contribution. Again, there are no contributions from the tangential 
and no- slip BCs. 

The relevant part of w^ is 

x^^ x"^ 19x1° 47a;9 x^ 2x'^ Ax^ 



4-'^=!::^^-^;::^ + :::^-:.:^ + ^ + 



95040 7920 28350 22680 420 315 135 



2x^ 58x^ 232a;3 12358a; \ , 
+ 7T7F - ^T7^ + 0110. R'n. 



45 ' 945 945 ' 31185 ' " ''^^^ 



10 „9 on^8 Q1^7 7„6 0^5 



X X 



29xS 31x^ 7x^ 2x^ 2x^ Ax'^ 344x\ , 



2700 270 1680 630 90 75 15 15 945 



4x^° 8x9 ^98 g 7 6 ^5 ^4 2x=^ 107x\ „2 .^2 

+ + + — + R^ cot ev^i 

14175 2835 1680 315 45 25 18 45 810 / 



Using the expressions for V2 and W3 in the incompressibility condition ( [B23|) , we obtain 



(<4) _ / x^^ x^2 19x" 47a;i° 

"^ ~ ~\ 1235520 ~ 95040 ^ 311850 ~ 226800 



Q 

X X 



4x^ x^ 58x5 5g^4 6i792;2^ 



3780 1260 945 135 4725 945 31185 / ' 



+ 77T7T -T7777 + T7T7-7T7T-777T + 7T7r- TTTTTTT ^ COt OV^rizz- 



V 124740 11340 1512 336 126 90 450 90 11340^ 

(B37) 

Following our standard procedure [see Eqs. (|78|) and (^)], the kinematic condition at x = 1 
is written as 

m + {wo + wi)riz 

= uo^rj + ui + ui^rj + U2 + U3 (x = 1) (B38) 
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where we have taken into account Eq. (0). Also, terms wirjz and uixt] which could normally 
be discarded in comparison with ui, must be kept because Ui{x = 1) contains rj^x with the 
factor 6 which is allowed to be small. Note that we have not shown terms of type Vorjy, virjy, 
V2i]y and W2i]z as they make smaller contributions. Also, wnxV = at x = 1. Using the 
expressions for the velocities (|§), (|73|), (|2|), ([T^), (|B22|), and (|B37p in Eq. 



we obtain 



5 4 295 

V - Y^Rvz + —Rcotev^T] - ^R\ 



+ Ar]r]z + —RVzz - ^ cot OV'^r] 
t 15 3 






23 

15^ 



i? - 2 cot ( 



5 2 6 

ivvz)z - 77^ cot ev% + ^R%zz - T cot ev^T] 

14 7 5 



331 „^2 

H R^ Vz 

168 ' 



1241483 
8108100 



R% 



477523 
2494800 



R^ cot eW^T],, = 0. 



(B39) 



In this equation, we do not shown the linear terms with derivatives of order greater than 
four. Also, we have not shown the nonlinear terms with the derivatives of odd orders, of 
the type rjzrjyy, rjyrjzy, {rjrjzz)z, etc. Using a procedure identical to the one of Refs. p^,p9 



(see also Ref. |3y), it can be shown that the nonlinear terms with odd-order derivatives 
are dispersive in nature, (i.e. they do not give rise to any change in the amplitude). In 
comparison with the linear dispersive term V^r^^, the nonlinear dispersive terms are small. 
Similarly, nonlinear terms with even-order derivatives can be shown to be dissipative. Among 
the nonlinear dissipative terms, the terms involving W^are of the type WrjV^ri^ Wrjy'V'^riy, 
etc. These terms are smaller than the linear dissipative term ~ WV^r], cot Orjyy, and 6r]zz by 
the factor A. However, the nonlinear dissipative term R{ririz)z may be ~ 6t]zz when 5 -C /2. 

APPENDIX C: GENERAL ANALYSIS OF ITERATION PROCEDURE 

We want to analyze the general nth step of the iteration procedure. Since, as was pointed 
out in the main text, this is of interest mainly for the near-critical conditions, R ^ Re = 
(5/4) cot 6, where the flow is essentially 2D, Y ^ Z, for simplicity we will consider from the 
outset a 2D flow {dy = v = 0) with cot 9 = 0. 
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From the preceding iteration steps, the corrections are known up to (including) Un-i, 
Wn-i, and Pn-2- We substitute into the exact NS problem the expansion 



where Wn is an unknown correction and all the other members of the sum are known from 
the preceding iteration steps, and similar expansions for u and p. The exact equations and 
boundary conditions for p^i are 

n—2 n— 1 /n— 1 \ /n— 1 \ 

Rp(n-l)x = -Yl ^Pkx + XI [^kxx + Ukzz " RUkt " R{wn - '2)Ukz] - RI^uA [^Ukx] 
k=0 k=0 \k=l J \k=l J 



and 



/n-l \ /n-1 \ 



erms containing Wn or -u^ 



n—2 n—1 

Rpn-i = -RPN -Y.RPk-2(l-rjl + r]t + --- + {-Ifrj^J') ^ ^kx 

fc=0 fc=0 



(CI) 



n-\ 



2 (r^, - r]l + r]l ) ^ {ukz + Wkx) + 2 (r^, - 7]] + r]l 



k=0 



WNx 



n-l 



— 2 fr]^ ~ ''?2 + ^2 ^ ■ ■ ■) X! ""^^^ "I" [terms containing w„ or m„] {x = h), 

k=0 

and the corresponding equations for the approximate Pn-i are written as 



(C2) 



RP(n-l)x = U(n-l)xx + U(^n-2)zz " RU(n~2)t — RiuiN — '2)U(n-2)2 



R^k=0,..,n-3 UkU(^n-3-k)x + WkU(^n-3-k)z 



(C3) 



and 



^n n~2 „n— fc— 1 

Rpn-1 - -RPn ^ - L ^P. (n-k-l)\ 



[{n-k-l)/3] 

E 

fc=0 



2(-i)^^: 



ra-3p-l „_fc_3p_i 
p^2p V^ „,("-fe-3p-l) 'I 



2 /l^ ""fci: 
A;=0 



(n — k — 3p — 1] 



53 



lin^k~3)/3] 

+ E 



p=0 



2(-i)^r' E 



n-3p-3 ^n-fc-3p-3 

{n-A:-3p-3) ^ 



■u 



fcz 



A;=0 



(n — A; — 3p — 3)! 



[(n-fc-2)/3] 

+ E 

p=0 



n-3p-2 n-fc-3p-2 

2(_l)Pr72p+l y- ^(n-fe-3p-2) ^ 



A:=0 



(n-/c-3p-2)! 



[(„_fe_4)/3] 

E 

fe=0 



?i--3p-4 n-fe-3p-4 
^p^2p+2 Y^ . ,(»-fc"3p~4) '/ 

fc=o '^ (n-A;-3p-4)! 



2(-l)^r^r^ E 



+ 2 (j], - r]l + r]l ) Wnx {x = 1) 



(C4) 



where the notation [■ ■ ■] above the summation sign J2 indicates that only the integer part of 
the expression inside the brackets needs to be considered. 

After the solution Pn-i has been found, the exact equations 

n-l 

Wnxx = E [-"^kxx - Wkzz + R{wn - 2)Wkz + RWkt 
k=0 



+RUkWNx + RP 



kz\ 



(n-l \ /n-l \ /n-l \ /n-l 

E "fe E ^fc^ +R[^Wk] i^Wkz 

k=0 ) \k=0 / \k=0 / \k=0 



+ [terms containing w„, «„, or Pr, 



(C5) 



n-l 

Wnx = -WNx - E (^fc^ + "fe^ 
A;=0 



n-l 



+2 Y^ {wkzVz - UkxVz) [^-vl 

k=0 



-1 



+ [terms containing w„, m„, w„ or p„] {x = h), 



(C6) 



and 
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Wn = {x = 0) 



(C7) 



lead to the system for w„ 



Wnxx = -W^n-l)zz + R{wn - 2)w(„_i)^ + RW(n~l)t 



+RU(n-l)WNx + RP{n-l)z + -RSfc=o,..,n-2 UkWl^n~2-k) 



+WkW{n-2-k)z 



(C8) 



"^nx — "^ Nx 



„n+l n-1 / ^n-k ^n-k-l \ [{n-k~2)/Z] 



[n + 



k=0 



{n-k)\ "' {n-k-l)\^ 



p=0 



n—3p—2 

X E 

A;=0 



W 



(n-k-Zp-2) {n-k-3p-2)\ V 



kz 



— U 



kx 



n—k—3p—2 



{n — k — 3p — 2) 



{x = 1), (C9) 



and 



Wn = {x = 0). 



(CIO) 



Finally, after the solution Wn is known, the system 



n— 1 

Unx = -Yl ("fc^ +Wkz) - Wnz " W(^ri+l)z 
k=0 



(Cll) 



and 



Un = {x = 0) 



(C12) 



leads to the problem for m^ 



Unx U!nz 



(C13) 



and 



Un = {x = 0), 



(C14) 
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solution of which completes the nth iteration step. In view of the exact equations for p^i, 
nTn, and u^, it is easy to justify by mathematical induction the equations for the approximates 
Pn-i, Wn, and which are the principal objects of our consideration here. 

Since each of these problems is a linear (nonhomogeneous) ODE with linear boundary 
condition, its solution is the sum of contributions from each term of the RHS of the equation 
(which contribution is the solution of the ODE with all other RHS terms put to zero and 
with homogeneous BC) and of each term of the free-surface BC (the latter contribution is 
the solution of the homogeneous ODE with the BC truncated to that single term). With 
this, it can be seen, also inductively, that the order of magnitude of each term of Wn is the 
product of certain powers of our four independent small parameters — A, 1/L^, R/L, and 
R/T — with the sum of the exponents being (n + 1), the same for all terms of Un (we will 
say that the term has the "order of smallness" n + 1); and, moreover, all possible types of 
terms of order {n + 1) are present in general, except for the terms of the type (i.e. order 
of magnitude) A^~^^. Essentially the same orders of magnitude are found in u„, since, from 
( |U13D , it is the sum of contributions of each term of w„ (differentiated in z, which gives the 



additional factor 1/L in the order of magnitude, and integrated in x, which does not affect 
the order of magnitude at all). Similarly, Rpn-i is made of all possible terms of the order 
n times the factor 1/L (in this case, without any exceptions, since the terms of type A^/L 
are present). One can note that the first three terms of the RHS (|C8|) give contributions 
which are smaller than w„_i by the factors l/L"^, R/L, and R/T, respectively, while the BC 
term W(^n-i)xxV contributes terms which are smaller by the factor A (including the terms 
whose types contain A"', such as A^R/L). A similar role for Rpn-i is played by the BC term 

RPn-2V- 

With the above order-of-magnitude structure of u„ being known, it is easy to see that 
each term of it is much smaller than some term of Un-i, by a small parameter 1/L^, R/L, 
R/T, or A. [It may seem at first sight that terms of type, say. A"' {R/L) present a difficulty: 
there is no A"' in Un-i (we omit the common 1/L). However, A^{R/L) = AlA^^^^R/L)], 
and the type A"'~^{R/L) is present in «„_!.] It follows that 
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Un < Un-l 

Thus, the iterative corrections of u are monotonously decreasing. (The same can be seen for 
the sequences {wn} and {pn}-) 

All Un are substituted into the kinematic equation, and the time-derivatives are elimi- 
nated by iterating the EE (similar to sec. |V|). Then we have, in addition to the terms of 
the GEE (|8^ ) (which are of order two), also linear and nonlinear terms of order three and 
more (not counting the factor l/L). Of these, the nonlinear dissipative (i.e. with the total 
number of 2;-differentiations being even) terms are of a type containing A to the power at 
least two, and containing R/L to the power not less than one (since there are at least two 
differentiations). Thus, the Ojvf-type of each of these terms contains all the type- factors of 
the (nondegenerate) nonlinear dissipative term ex R{rf)zz, ~ iXl L)A^{R/ L), and, except for 
this principal nonlinear dissipative term, also some additional small factor, the product of 
powers of the small parameters A, 1/-^^, and R/L, which makes that term negligible. There- 
fore, we do not have to retain any nonlinear dissipative term other than the principal one 
(which is of order three). Similarly, of all the dispersive terms (both linear and nonlinear), 
only the (nondegenerate, second-order) term oc ri~,zz is retained. 

As we noted in the main text, all the linear dissipative terms could be similarly neglected 
as compared with the linear second-derivative one, but the latter is degenerate if 5 -C i?. 
Therefore, we should consider these linear terms. If the type of a linear term does not 
contain a factor which is a power of R/L, that term is just a (2n -|- l)-derivative one, clearly 
a dispersive term. Thus, the dissipative linear terms we are interested in must have the 
(/th-derivative) form R^ri^^\ which leads immediately to Eq. ( |101|) . 
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FIGURES 
FIG. 1. Evolution of energy illustrating that the solutions of Eq. (p6|) remain bounded. 



FIG. 2. Evolution of energy indicates that the solutions of Eq. (|104|) blow up. 
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